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We study a QCD sum rule analysis for an exotic tetraquark ud bar[s] bar[s] of JP=0+ and I 
VQ \ — 1. We construct q q bar[q] bar[q] currents in a local product form and find that there are five 

f*^ . independent currents for this channel. Due to high dimensional nature of the current, it is not easy 

to form a good sum rule when using a single current. This means that we do not find any sum 
£"\j | rule window to extract reliable results, due to the insufficient convergence of the OPE and to the 

exceptional important role of QCD continuum. Then we examine sum rules by using currents of 

linear combinations of two currents among the independent ones. We find two reasonable cases that 

predict a mass of the tetraquark around f .5 GeV. 
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The history of exotic hadrons is rather long. But the recent experimental observations have triggered tremendous 

amount of research activities Ju 0, El ■ Among them the report on the pentaquark + from the LEPS group in 2002 

was the most influential one |^ , partly because + is a genuine exotic state of the quark content uudds. It also has 

unusual properties such as a light mass and a very narrow width. Its existence is, however, now questioned, which 

i-G ' should be confirmed in the future experiments [5j- 

*~H. Turning to mesons, though not genuine exotic states, -XY3872) and D s (2317) are found to have properties that seem 

O ,' difficult to be explained by a conventional picture of qq [3, 0, M, lj EH HH 03 Rather, they could be considered to 

have a significant amount of multi-quark components. Historically, tetraquark mesons were investigated long ago as 

rG an attempt to explain relatively light masses and excess of states in scalar channels [OL HH Ha Ho . Il7| . Just as in 

^. , the exotic baryons, it is interesting to consider genuine exotic states in the meson sector whose minimal component 

is qqqq. Tetraquark states of udss component have been studied as candidates of such exotic states. Since they may 

be obtained by replacing one of ud diquarks in + by an s a ntiq uark, similarities between + and udss have been 

discussed, though precise analogy is a dynamical question [lj, Oil 123 • 

In the former studies, the tetraquark udss of J p = 1 + was investigated in detail, where it was shown that the state 
has a relatively low mass and a narrow width decaying into K*K in the flux tube model |22j. The narrow decay width 
is associated with the fact that KK channel is forbidden due to the conservation of parity and angular momentum, 
which partly motivated the study of the 1 + channel. 

In principle, it is also possible to study other channels of the udss tetraquarks |2lL l22l |23| . From a naive point of 
view of mass, it is natural to investigate + scalar states. In contrast to qq mesons, the tetraquark does not need 
orbital excitation to form the quantum number + , but all quarks may occupy the lowest state. In this case, it is 
shown that the tetraquark should have isospin one 1=1. This is the object that we would like to study in this paper. 
We perform QCD sum rule analyses for the scalar (J p = + ) and isovector (/ = 1) exotic tetraquark udss. We 
attempt a rather comprehensive analysis in which we will pay special attention to the structure of the interpolating 
fields (currents). First, we find that there are five independent interpolating fields for the tetraquark. We show this 
by constructing the tetraquark currents in terms of diquark fields {{qq){qq)) and mesonic fields ((qq) (qq)), where qq 
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can be both color singlet and octet. We then consider two-point correlation functions first by using a single current 
of various types. It turns out that many of them do not achieve a good sum rule. Therefore, we attempt linear 
combinations of two independent currents. This method was first proposed in Ref. [2J|. We then find that there are 
several cases with good Borel stability, indicating the mass of the tetraquark around 1.5 GeV. We also investigate the 
reliability of the sum rule not only from the Borel stability but also from the dependence on the threshold value and 
the amount of the pole contribution in the total sum rule. We also mention the convergence of OPE. 

The difficulties to make a good sum rule for exotic particles of high dimensional operators were nicely discussed in 
a recent work by Kojo et al. |25J. They proposed a sum rule using a linear combination of two-point functions rather 
than currents in order, for instance, to suppress large contributions from low dimensional terms that are irrelevant to 
non-perturbative properties of hadrons. They have successfully achieved a good sum rule that satisfy the necessary 
requirements. In our present study, our strategy is different from theirs, but the consideration along their idea is 
certainly important in the discussion of the tetraquark also. 

This paper is organized as follows. In section 2, we establish five independent currents in diquark-antidiquark and 
meson-meson (actually meson-like) constructions. Some relations among various currents will be discussed. Section 

3 is the main part of this paper, where we perform sum rule analyses using various tetraquark currents constructed 
in section 2. We study the sum rule of a single current and then consider linear combinations of currents. Section 

4 is devoted to summary. In Appendix, we discuss the equivalence and relations between the currents of diquark- 
antidiquark and meson-meson constructions. 

II. INDEPENDENT CURRENTS 

Let us consider currents for the tetraquark udss having J p = + . Here we consider only local currents. To write 
a current, Lorentz and color indices are contracted with suitable coefficients (L^^) to provide necessary quantum 
numbers, 

V — lj nupcr s a s b u c a d ! \ L ) 

where the sum over repeated indices (fi, !/,■■■ for Dirac spinor indices, and a, b, ■ ■ ■ for color indices) is taken. 

For the Dirac spinor space, using possible diquark and antidiquark bilinears [261 I27L l28ll29| . there are five indepen- 
dent terms 

Sabcd = {SalbCsl){u^C^d d ) , 
Vabcd = (s a "ff_i"f5Csl)(u^Cj tJ '"f5dd) , 

T a bcd = i.s a a^Csl){u T c Ca^d d ) , (2) 

Aabcd = (s Q 7 /1 Cs^)(uJC , 7 M dd) , 

Pabcd = {SaCsl){u T c Cd d ). 

Here, color indices are not yet specified. For the diquark and antidiquark pair, color structures providing a color-singlet 
tetraquark are 3 (g> 3 and 6 <g> 6, which we will denote by labels 3 and 6 for short. 
Therefore, we have altogether ten terms of products 

{S®V®T®A® P} Lorentz <S> {3 © 6} C olor ■ (3) 

However, half of them drop due to the Pauli principle. For instance 

-P3 = Phorentz ® ^Color (4) 

= e a bc(s b Cs^)e a b' c '(Ub'Cd c ') = 0. 

Eventually, we end up with five independent currents 

^6 = (s Q 7 5 Csf )(m^ C754) , 

V 6 = (salpl5Csl)(ulcri 5 d b ), 

T 3 = {s a ap V Csl){u T a Ca^db) , (5) 

A z = {s al ^Csl){u T a Crd b ), 
P & = {s a Csl){u T a Cd b ). 



In the non-relativistic language, these five terms correspond to combinations of diquarks and antidiquarks 

[( 1 5o)( 1 ^o)] 0+ , [( 3 Si)( 3 Si)] 0+ , [( 1 Fi)( 1 Pi)] 0+ , [( 3 Po)( 3 Po)] 0+ , [( 3 A)( 3 Pi)] 0+ ■ (6) 

Another possible piece of Z P^ is irrelevant, since the five bi-linear forms q T Tq (T = S, V, T, A, P) can only have spin 
j < 1, while the 3 P2 diquark has j = 2. 

Finally we consider the flavor structure. The ss antidiquark is symmetric in flavor, and hence belongs to the 
symmetric representation 6/. If the other ud diquark belongs to 3/, and so isospin I = 0, the diquark and antidiquark 
will have different flavor symmetry. But they should have the same color and spin symmetries for composing a color- 
singlet scalar tetraquark. Considering the Pauli principle, they must have different parity, and hence their combination 
is a negative-parity scalar tetraquark. Accordingly, the other ud diquark also belongs to 6/, and so isospin 1=1. 
Among the irreducible representations of the tetraquark 

6<g>6 = l©8©27, 

S = +2 and 1=1 states arc in the 27 representation of SU(3)f, which is the flavor structure of the present tetraquark. 
As shown in Fig.^ three isovector states of the 27/ are uuss, l/\/2(ud + du)ss and ddss. 



dd s s 




5 = 2 



FIG. 1: SU(3) weight diagram for 27, where the locations of three tetraquark components of S = 2 and 7=1 are shown. 

We have constructed five independent currents using diquark and antidiquark combination. We refer to this as 
the diquark construction. Similarly, we can also construct the tetraquark currents using qq combination (mesonic 
construction). Obviously, there are ten combinations of the Dirac (5, V, T, A and P) and color (1 and 8) spaces: 



5*1 = (s a u a )(s b d b ) , 

Vi = {s a ~f^u a )(s b -f^d b ) , 

T\ = {SaCr,j,uUa){Sb^ V db) , 
M = (Sa7/*75W )(S67 M 7 5 4) 

Pi = (s a j5U a )(s b J 5 d b ) , 



Ss = {s a \ ] b u b ){s c \^ d d d ) , 

Vs - {s a -i^\ n ab u b )(s cl ^\ n cd d d ) , 

T 8 = (s a a^y: b u b )(s c <j^X n cd d d ) , 

A 8 = {Salnl<iKb U b){Sc^l<iKddd) , 
Ps, = (Sa75A" b M&)(s c 7 5 A" d d d ) , 



(7) 



where subscripts 1 and 8 denote color singlet and octet representations, respectively. Unlike the diquark construction, 
all the ten currents in Eq. Q remain finite. However, it is possible to show only five of them (in fact any five of them) 
are independent. The proof of this and various relations among different currents are discussed in Appendix. 1X1 



III. QCD SUM RULES ANALYSIS 
A. Formulae of QCD Sum Rule 



For the past decades QCD sum rule has proven to be a very powerful and successful non-perturbative method |3J . l31| . 
In sum rule analyses, we consider two-point correlation functions: 



n( g 2 ) = i f d A xe iqx {U\Tr]{x)rf{<S)\0), 



(8) 



where r\ is an interpolating current for the tetraquark. We compute H(q 2 ) in the operator product expansion (OPE) 
of QCD up to certain order in the expansion, which is then matched with a hadronic parametrization to extract 



information of hadron properties. At the hadron level, we express the correlation function in the form of the dispersion 
relation with a spectral function: 

n(p)= r p{ ° ] . da, (9) 

Jo s-P - *£ 

where 

p(s) = ^^--^nXOMnXnl^lO) 

II 

= fx s ( s - M x) + higher states . (10) 

For the second equation, as usual, we adopt a parametrization of one pole dominance for the ground state X and a 
continuum contribution. The sum rule analysis is then performed after the Borel transformation of the two expressions 
of the correlation function, © and @ 

/>oo 

rt all \M%) = B M% W{p 2 ) = J e- s ' M Bp{ s )d S . (11) 

Assuming the contribution from the continuum states can be approximated well by the spectral density of OPE above 
a threshold value so (duality), we arrive at the sum rule equation 

U(M B ) = f xe -M 2 x /Ml = r° e s/Ml p ^ ds (12) 

Jo 

Differentiating Eq. (|12[1 with respect to -rjT an d dividing it by Eq. (|12|l . finally we obtain 

M B 

2 _ J So e- s / M2 Bsp(s)d S 

X ~ Ce-^lp(s)ds • (13) 

In the following, we study both Eqs. (|12J) and (|13|l as functions of the parameters such as the Borel mass Mb and the 
threshold value so for various combinations of the tetraquark currents. 

B. Analysis of Single Diquark Currents 

In this subsection, we perform a QCD sum rule analysis using the five diquark currents, Sq, Vq, T 3 , A 3 and P§, 
separately. Let us first outline briefly how we performed the OPE calculation. For illustration, let us take Pq. Then 

n(g 2 ) = ifd^xe^iOlTPe^P^OM 

= Tr[C{Sf {x)) T CSf {x)]Tr[S:' a {-x)C{S b : h (-x)) T C] (14) 

+Tr[C{Sf (x)) T CSf (x)]Tr[S»: a {-x)C{S«'\-x)) T C]. 

For the quark propagator, we use 

iSf(x) = (0\T[ q a (x)q b (0)}\0) 

AXab : \n i sab A ah T 2 

= 2 ^^^^^^^ + ^)-T 2 -^ + T92-^^ ) (15) 

5 ab m q iS ab m„ , . iS ab m 2 a 

— , n n H 7^ {qq)X+ a 9 9 X ' 

Att 2 x 2 48 8n 2 x 2 

The two-point function is then divided into three parts: 

1. Terms proportional to S ab (a, b being color indices), where no soft gluon is emitted. The lowest term of this kind 
is the continuum term. 

2. Terms containing one A a h (color matrix), where one soft gluon is emitted. The lowest terms of this type contain 
condensates such as (gqaGq) (q — u and d) and (gsaGs). 



3. Terms containing two A a b's, where two soft gluons are emitted. The lowest terms of this type contain the 
condensate (g 2 G 2 ). 

We have performed the OPE calculation for the spectral function up to dimension eight, which is up to the constant 
(s°) term of p(s). Actual computation is very complicated. We have performed this calculation using Mathematica 
with FeynCalc |32(. Mathematica programs are available from the authors. The results are 



a 4 m s 2 s 3 t m£_ _ m s (ss) _ (g 2 GG) \ 

Pse( s ) : = 614407r 6 3072tt 6 V256tt 6 192tt 4 12288tt 6 / 



/ m 3 (ss) m 2 (g 2 GG) m s {gsaGs) (qq) 2 (ss) 2 \ 

\ 32tt 4 4096tt 6 64tt 4 24tt 2 24tt 2 / S l ' 

m 2 (qq) 2 m 2 (ss) 2 (gq)(gqaGq) m s (g 2 GG)(ss) (ss)(gsaGs) mj(g 2 GG) 
Yin 2 48tt 2 24tt 2 1536tt 4 24tt 2 2048tt 6 

„4 K™ 2 3 , m 4 ™ I-e \ z/„2. 



5m s s 3 f m s m s (ss) 5(g GG)\ 
PV ^ S > = 15360tt 6 ~ 1536tt 6 + V64^ + 24tt 4 + 6144tt 6 ) S 



.2 

mj}<5s) llm 2 (. 9 2 GG) m s (. 9 scrGs) (gg) 2 



8tt 4 2048tt 6 32tt 4 12tt 2 12tt 2 . 

2m 2 s (qq) 2 m 2 (ss) 2 (qq)(gqcrGq) 7m s (g 2 GG)(ss) (ss)(gsaGs) 



(17) 



3tt 2 12tt 2 12tt 2 768tt 4 12tt 2 

n _ s 4 rn s 2 s 3 /3m s 4 m s (ss) (g 2 GG) x a / 3mf(£g) m 2 (.g 2 GG) \ 

PT3lSj 5120tt 6 128vr 6 V647T 6 16tt 4 1536^ 6 / V 8tt 4 256vr 6 / 

m 2 (gg) 2 m 2 (ss) 2 m,(g 2 GG)(ss) mj(g 2 GG) 

it 2 4tt 2 192tt 4 256tt 6 ' l } 

A m 2„3 , ^ 4 /„2/ 

„2 



= S 4 m s V /_m/_ (g 2 GG) x 

pj43lSj 30720tt 6 1024tt 6 Vi28tt 6 61447T 6 /' 



m^(^) 3m 2 (ff 2 GG) m s (gsaG S ) (g'g) 2 (ss) 2 y 
16tt 4 2048tt 6 64tt 4 24tt 2 24tt 2 / 1 

[ss) 2 (qq)(gqaGq) m s (g 2 GG)(ss) (ss)(gsaGs) 



(19) 



24tt 2 24tt 2 256tt 4 24tt 2 

^ *» 2 - 3 ' - 4 m s (ss) (g 2 GG)\ 



s m s s / rr 

PP6(S) = 61440tt 6 ~ 1024tt 6 + V25 



l.s 2 



6tt 6 64tt 4 12288tt 6 /' 
/ m 3 (ss) 3m 2 s (g 2 GG) m s (gsaGs) (gg) 2 (ss) 2 \ , . 

V 32tt 4 4096tt 6 64tt 4 24tt 2 24tt 2 / S l ' 

m 2 (qq) 2 rn 2 (ss) 2 _ (qq)(gqaGq) _ m s (g 2 GG)(ss) _ (ss)(gsaGs) _ mj(g 2 GG) 
4tt 2 48tt 2 24tt 2 512tt 4 24tt 2 2048tt 6 

In these equations, g represents auord quark, and s represents an s quark, (qq) and (ss) are dimension D = 3 quark 
condensates; (g 2 GG) is a D = 4 gluon condensate; (gqaGq) and (gsaGs) are D — 5 mixed condensates. From these 
expressions, we observe the followings: 

• The coefficients of the lowest dimension, or of the leading term in powers of s, have the relations c^ 6 = c P ^ and 
c^g = 1/2c V q. These are the consequences of chiral symmetry at the perturbative level |33l |. 

• As empirically known, the terms of quark condensates have important contributions to the sum rule. 

For numerical calculations, we use the following values of condensates paL l3o . 1371 J38L l39l PlOl l41| : 

(qq) = -(0.240 GeV) 3 , 

(ss) = -(0.8 ± 0.1) x (0.240 GeV) 3 , 

(g 2 s GG) = (0.48 ± 0.14) GeV 4 , 

m s (2 GeV) = 0.11 GeV, (21) 

(g s qaGq) = -M 2 x (qq) , 

M 2 = (0.8 ± 0.2) GeV 2 . 

As usual we assume the vacuum saturation for higher dimensional operators such as (0|gggg|0) ~ (0|gg|0)(0|gg|0). 
There is a minus sign in the definition of the mixed condensate (g s q<jGq), which is different with some other QCD 
sum rule calculation. This is just because the definition of coupling constant g s is different |34l l35| . 
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FIG. 2: Borel transformed correlation functions IIs6, 11^6, IIt3, IIa3 and IIp6 as functions of Borel mass square, in units of 
GeV 10 , for threshold value s = 3 GeV 2 . 



In Fig. we show all five Borel transformed correlation functions Ii{M%) (the LHS of the Eq. I[12)l) as functions of 
Borel mass square for threshold value so = 3 GeV 2 . From the definition of I|1U|) . the LHS should be positive definite 
quantities. In practical calculations, however, the positivity may not be necessarily realized, if the OPE up to finite 
terms does not work due to insufficient convergence of the OPE. In the present analysis, we find that among the 
five cases, two functions of Vq and Pq currents show such a bad behavior. Therefore, the QCD sum rules for these 
two currents are not physically acceptable. The correlation functions of A3 and Sq change the sign from negative to 
positive values. But the sum rule values take positive values for M% ~ several GeV 2 . 

The tetraquark currents Sq and A3 are constructed by diquark fields which correspond to 1 Sq and 3 Si in the 
non-relativistic language, where the two quarks can be in the ground state s-orbit. In contrast, the currents Vq and 
Pq correspond to linear combinations of 3 Pi, and 3 Po, respectively, where one of the two quarks is in an excited 
p-orbit. The T3 current is a linear combination of 3 S\ and 1 P\. Therefore, we verify an empirical fact that the 
sum rule constructed by currents having the s-wave components in the non-relativistic limit works better than those 
dominated by p-wave components. For completeness, we show the LHS with numerical coefficients for the three better 



cases: A3, T 3 and Sq 



n 



(all) 
A3 



8.2 x 10- 'M 



-7»rl0 



7.4 x l(T 8 Ml 



1.6 x 10~ 7 M% + 1.8 x 10~ 6 M% 



1.1 x KT 6 Ml 



n^ Q 3 = 4.8 x 10~ 6 Af^° - 5.9 x KT 7 Af| - 9.1 x lCT 7 Af| 



3.4 x 10~ 8 Af| + 2.4 x lCT 7 Af| 



n 



(all) 
S6 



4.1 x Kr 7 AfJ, 



2.5 x 10 _8 Ml 



5.1 x 10 _8 M| + 1.8 x 10~ 6 Ait - 1.1 x 10~ 6 Mp 



(22) 



From these expressions, we observe that the convergence of the current T 3 seems better, while the convergence of the 
currents A3 and Sq is not very good in the region 1 < AfJ, < 2 GeV 2 . They can only converge at A/| ~ 3 GeV 2 . 

To determine the mass, we need to fix the two parameters: the threshold value sq and the Borel mass square Mg. 
For a good sum rule, the predicted masses should not depend on these two parameters strongly with sizable pole 
contribution (Borel window) . In Fig. |21 we show the masses of the tetraquark as functions of the Borel mass for 
several threshold values sq (Borel curves). We observe that the Borel mass dependence is somewhat strong for the 
currents S& and A3 in the region 1 < Af| < 2 GeV 2 , which is expected to be a reasonable choice of the Borel mass. 
For these currents Sq and A3, however, we see that the minimum occurs at around 3 GeV 2 when so is varied in the 
region M\ > 1.5 GeV 2 . (For the current Sg, the mass of Sq = 2 GeV 2 is far above the region shown in the figure.) 
For this reason, we consider that so = 3 GeV 2 is a reasonable choice which we will mainly use for the estimation of 
the mass of the tetraquark in the following sum rule analyses. At this so value, the mass of the tetraquark turns 
out to be about 1.6 GeV. For the T 3 current, the Borel stability seems better. The result, however, depends on the 
threshold value so to some extent. However, it is interesting to see that the mass of the tetraquark is about 1.6 GeV 
when s - 3 GeV 2 . 
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FIG. 3: Mass of the tetraquark calculated by the three currents Se, A3 and T3, as functions of the Borel mass square A/|, for 
several threshold values so = 2, 3, 4 and 6 GeV 2 . 



To see the amount of the pole contribution, we define the quantity 



Pole contribution = 



r 

j 4 



,-s/Mf 



p(s)ds 



f- 



-a/Ml 



p(s)ds 



(23) 



As shown in Table 1, the pole contribution of the diquark currents A3, T 3 and Sq are not very large; at Mg ~ 1 
GeV 2 they are of order 10 %. This is a general problem of the QCD sum rule when multi-quark currents are used. 
Therefore the results so far might be doubtful. 



TABLE I: Pole contributions of various currents. The threshold value so = 3 GeV is used. 



Diquark Current 
A 3 T 3 S 6 



Mesonic Current 
V a T 8 



Mixed Current 

Vi V'2 



0.7 GeV 

1 GeV 2 

2 GeV 2 



0.17 0.11 0.10 
0.04 0.01 0.05 



0.54 
0.09 



0.23 
0.02 



0.60 
0.30 
0.03 



0.49 
0.22 
0.02 



From the analysis of the single current of the diquark construction, we expect that the mass of the tetraquark is 
about 1.6 GeV, although the stability against the variation of both the Borel mass and the threshold parameter is 
not simultaneously achieved. Furthermore, the pole contribution is rather small. As we will see, however, a suitable 
linear combination will improve them. 



C. Analysis of Single Mesonic Currents 



In this subsection, we perform QCD sum rule analysis using the ten mesonic currents, Si,s 
Pi, 8) separately. Here we only show two important spectral densities: 



Vi,s, Ti «, Ai,s and 
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54tt 2 
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864tt 2 
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+ 
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m 3 (gqGq) 

1447T 4 
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432tt 2 



(24) 



'»: 



18432tt 6 
/ 5m 3 - 



( 



2304tt 6 V3847T 6 
ss) 31m 2 (g 2 GG) 



+ 



1728tt 2 
5m s (ss' 



18tt 2 

<{g 2 GG) 

9216tt 6 



(ggGg) (I 
864tt 2 



31 (ff 2 GGh 2 
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5m 2 (qq) 2 



5m 2 (ss) 2 
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9216tt 6 



18tt 2 



72tt 2 



(25) 



As shown in Fig. 0J we find that among the ten correlation functions, only two correlation functions for the currents 
Vs and Ts show good behavior with having positive values. 

The currents V\ , Vs , Pi and P$ are constructed by mesonic fields (either color singlet or color octet) which correspond 
to 3 S\ and 1 Sq in the non-relativistic language, where two quark-antiquark pairs can be in the ground state s-orbit. 
Their spectral densities then show similar behavior to 5"6 and A3 in the previous subsection. In contrast, Si, Ss, A± 
and Ag correspond to linear combinations of 3 Pq and 3 P\, respectively; T\ and Ts currents are the conrbinations of 
3 Sr and l P x . 

From the above argument, we might expect that six currents, Vi, Vs, -Pi, -Ps: T\ and Ts would work. However, we 
found that the Borel transformed correlation functions calculated by the currents Vi, Pi, Ps and Ti take negative 
values and therefore, they must be abandoned. Now there remain only two better currents Vg and T 8 in the mesonic 
construction. This is the reason why we have shown their spectral densities in 124fl and (|25|) . Using the numerical 
values of various condensates (|21|l . we find the Borel transformed correlation functions 



(all) 

V8 

(all) 

T8 



2.3 x l(T 7 Af|, - 2.6 x KT 8 ATl + 9.1 x l(T 8 Af|, + 3.5 x 10" 7 Af| 



1.4 x 1Q-°M 



6JM-10 



1.7 x lQ- 7 Ml 



1.2 x IQ- 7 M% 



4.3 x 10 



" 9 M 4 



4.9 x 10 
4.9 x 10 



" 8 Mi 
" 8 Mi 



(26) 



From these equations, we find that better convergence is achieved for Ts than for Vs in the region 1 < M| < 2 GeV 2 . 
The pole contributions are significantly improved as shown in Table 1. 
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FIG. 4: Borel transformed correlation functions IIsi, flss, HVi, Ilvs, Uri, IIts, Hai, IIas, IIpi and IIp8 as functions of Borel 
mass square, in units of GeV 10 , for threshold value so = 3 GeV 2 . 
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In Fig. |SJ we show the masses of the tetraquark currents V$ and Tg as functions of the Borel mass for several 
threshold values s (Borel curves). As in the case of T 3 current, the Borel stability seems good but the result depends 
on the threshold value sq. However, once again, if we take the threshold value at so ~ 3 GeV 2 , the mass of the 
tetraquark turns out to be reasonable, though the precise values are slightly smaller: the mass of Tg ~ 1.5 GeV and 
the mass of Vk ~ 1.4 GeV. 
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FIG. 5: Mass of the tetraquark calculated by the currents Vg(Left) and Ts (Right), as functions of the Borel mass square M B , 
for several threshold values so = 2, 3, 4 and 6 GeV 2 . 



D. Analysis of Mixed Currents 

In order to improve the sum rule, we attempt to make linear combinations of independent currents for both diquark 
and mesonic currents. Since linear combinations of five currents contain ten mixing angles, the full consideration with 
these ten parameters is rather cumbersome. Instead, we make a linear combination of two currents J\ and Ji (any 
two from the independent currents), r\ = cos6>Ji + sin6>J2, where 9 is a mixing angle. Then the correlation functions 
are written as 

(r)rf) = cos 2 d{JiJ\) + sin 2 6»(J 2 4) + cos 9 sin 9(Ji J.J) + cos sin 0<J 2 jJ) . (27) 

The mixing is chosen with the following requirements: 

1. The OPE has a good convergence as going to terms of higher dimensional operators. 

2. The spectral density becomes positive for all (or almost all) s values, and then n(Af|) becomes positive for all 
Borel mass and threshold values. 

3. Pole contribution is sufficiently large. 

We have tried various combinations of two currents to realize good sum rules. While doing so, we have realized 
that the diquark currents are more independent than the mesonic currents. This means that the cross terms of H27fl 
have only a minor contribution for diquark currents, while they have a large contribution for mesonic currents. 

According to the requirement (1), we would like to make a linear combination such that the highest dimensional 
(eight) term is suppressed. For diquark currents, we find it convenient to take two combinations: 



rj = cos 9 A3 + sin 9V§ , 
f = cos 6P & + sin 9S 6 . 



(28) 
(29) 



By choosing cot 9 ~ y/2, we find that the term of dimension eight of (|28|l is suppressed, while for cot 9 ~ 1, the term 
of dimension eight of (|29|l is suppressed. The Borel transformed correlation function of (|29|l n^(Af^), however, takes 
negative values. Therefore, this current should be rejected for the sum rule analysis. In this way we are led to the 
current r) of J3SJ). From now on, we will denote r\ — > x\\. 

For the mesonic case, it turns out that the cross term contributions are large. Accordingly, we attempt a complex 
angle to improve the sum rule analysis. By choosing t\ — 0.91, £2 = —0.41, we construct a current: 



<Yl 



5i + (ti + i* 2 )Pi 



(30) 



The numerical Borel transformed correlation functions are 



(all) 

1 

(all) 



1.1 x 10 _6 Afi° - 1.3 x 10 _7 Afl + 4.8 x 10 _7 A/£ - 2.0 x 10~ 8 Af£ + 5.2 x 10~ 9 Ml 



B ■ 



5.0 x 10'' Mg - 6.0 x 10" S M| + 8.4 x lO^M^ - 2.2 x 10~ S M^ + 8.3 x 10~ M M|, . 
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(31) 



which may be compared with the previous results of l|22|l and (|26(l . It looks that the convergence of the series is 
improved significantly. Therefore, we can choose a smaller Borel mass square down to Af§ > 0.7 GeV 2 , where the 
pole contribution will be further increased up to around 50 %, and the convergence is still maintained. 
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FIG. 6: Mass of the tetraquark calculated by the mixed currents ?7i(Left) and 772 (Right), as functions of the Borel mass square 
M% for several threshold values so = 2, 3, 4 and 6 GeV 2 . 



In Fig. H3 we show the mass calculated from 771 and 772 as functions of the Borel mass square for several threshold 
values So ■ The Borel stability is improved from the cases of the single currents. From these figures, we might think 
that there is still a substantial sq dependence. However, this dependence will be largely reduced if we choose a small 
Borel mass, where the pole contribution is sufficiently large. In Fig. [7] we show the mass calculated from 771 and 772 as 
functions of the threshold value for several Borel masses. When Af| = 0.7 GeV 2 , the curve is very stable. Moreover, 
the pole contribution is around 50 %, and the convergence is still maintained. Therefore, we obtain a very good sum 
rule, where we find the mass calculated from the two currents 771 and 772 is about 1.5 GeV. As the Borel mass increases, 
the pole contribution decreases, and accordingly, the threshold dependence becomes bigger. 
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FIG. 7: Mass of the tetraquark calculated by the mixed currents 771 (Left) and 772 (Right), as functions of the threshold value so 
for several Borel mass square Mg — 0.6, 0.7, 1 and 2 GeV 2 . 



Finally, in order to summarize our analysis, we show in Fig. [S] masses of the tetraquark calculated by several 
reasonable currents used in the present study as functions of the Borel mass square at s = 3 GeV 2 . They are 5*6, A 3 
and T 3 for the diquark construction, T 8 and V% for the mesonic construction, and rji and 7/ 2 for the mixed currents. 
The plots are extended to a wider region of M^ up to 4 GeV 2 , where the masses predicted by different currents tend 
to a same value. We verify once again a good Borel mass stability for the mixed currents, while some of the single 
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currents show good stability also (T3,T§ and Vs). The mass values varies slightly, while we expect the mass of the 
tetraquark around 1.5 GeV. 
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FIG. 8: Mass of the tetraquark calculated by the currents jji, r\i, A$, Se, T3, Vs and Ig, as functions of the Borel mass square 
Ml in the region 0.7 < M% < 4 GeV 2 , for threshold value s = 3 GeV 2 . 



IV. SUMMARY 



We have presented a QCD sum rule study of the udss tetraquark of J p = + and 1=1, both in the diquark 
((qq)(qq)) and mesonic ((qq)(qq)) constructions. We have found that in this channel of tetraquark, there are five 
independent currents, which is shown both in the diquark and mesonic constructions. For each single current, we 
have tested the sum rule analysis, but it is found that not all of them provide a good stability. 

As an attempt to improve the stability of the sum rule, we have considered linear combinations of independent 
currents. In order to simplify the analysis, we took a superposition of various combinations of two currents. Among 
them, we have found two cases that lead to good sum rules, where we investigated Sq (threshold value) and Mb (Borel 
mass) dependence, and convergence of OPE. A good Borel stability is achieved in the region 0.7 < M B < 4 GeV 2 . 
In order to obtain a large enough pole contribution (50 %) and reduce the threshold value dependence, we have to 
reduce the Borel mass. However, to maintain the convergence of OPE, we can not reduce it too largely. When Borel 
mass square M B is around 0.7 GeV 2 , we get a very good QCD sum rule, where the mass of the tetraquark turns out 
to be around 1.5 GeV. 

Despite the seemingly good Borel mass stability, we think that we should investigate the following points more 
carefully. For instance, estimation of higher dimensional terms of 0(1/ s) could be important. Although we are able 
to construct the two mixed currents such that the higher order contributions (in the present calculation of OPE) of 
dimension six and eight terms are suppressed, the question still remains concerning even higher order contributions. 
Another question is the contribution of KK scattering states, since the mass of the tetraquark is around 1.5 GeV, 
and it can fall apart into the KK states. Such a contribution can be estimated by using the method proposed in 
Refs. J42,|4Jj. These will be further investigated in the future work. 
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APPENDIX A: FIVE INDEPENDENT CURRENTS IN (qq){qq) BASIS 
We attempt to write a diquark current of {J5J as a sum of (qq) mesonic pairs (q = u, d, s), 

L^pa X S%S%U p c dZ = J2 ClC%(8 b \ n r i dd)(8 a \ n T j Uc) , 



(Al) 



where I\ are the five Dirac matrices and X n (n — 1, ■ • • ,8) are color matrices forming color singlet and octet states 
out of 3 x 3. Therefore, in IjAljI . the sum runs over ten terms of five I\ matrices and two A™ combinations. They are 



Si = (s a u a )(s h db) , 

Vi = (8 a ,j li u a ,)(sb'y'*db) , 

P\ = {s a ^^u a ){s b cr^d b ) , 

M = ( S a 7m 75«a)(Sf> 7^754) , 
Pi = {SalbUa){Sblbdb) , 

where in the octet representation inner product of A™ [n = 1, 
currents have different properties: 

Si:(J F =0+ 8/,l c ), 
Vi:(J p = l-,8/,l c ), 
T 1 :(J p = l+&l-,8 / ,l c ), 
A 1 :(J p = l+,8 f ,l c ), 
P i: (J p = 0-,8/,l c ), 



S 8 = (s Q A" & u b )(s c A" d d d ), 

V S = {Sa1nKbUb){Scl' M Kddd) , 

T 8 = (s a a^X: b u b )(s c *^X> d d d ), 

A = {SalnlsKb U b)(ScYlsKddd) 
P S = {Sal^Kb U b){s cl5 \ n cd d d ) , 



(A2) 



S's 
V 8 

As 

Pr 



' is taken. The quark-antiquark pairs in different 

(J p = 0+,8 /; 8 c ), 
(J p = l-,8 /; 8 c ), 
(J p = l+&l-,8 / ,8 c ), 
(J p = l+,8 / ,8 c ), 
(J p = 0-,8 / ,8 c ). 



In order to establish the five independent currents, first we change their color structures 

{s a u b )(s b d a ) = -(s a u a )(s b d b ) + -(s a Ub)(s c dd)X a bKd , 
(s a Ud)(s c db)X a b\cd = -£-(s a u a )(sbdb) - ^{s a Ub){s c dd)Xab\cd ■ 
Then we use the Fierz transformation J44J 

-z(s a u a )(sbdb) + ■^(s a Ub){s c d d )X ab Xcd 



(s a u b )(s b d a ) 

-{(SaU a )0 

(s<j7m 75 Wa) (^7^754) + (s a 75"a)(sfc75<4)} 



j{(s a Ua)(s b db) + (Sa7M U a)( S fc7' i 4) + ~ (SaV U «) (S6«7 / "'dj ) ) 



(A3) 



(A4) 
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We obtain 10 equations in all 

-Si + ^s 8 = -i{5 x + v, + ir a - a x + a} , 

y 5i - i5 8 = -i{S 8 + y 8 + ip 8 - A 8 + P 8 } , 
^Vi + iy 8 = - ~{4Si - 2Ki - 2i4i - 4PJ , 

Y v i - l v * = -\i AS * - 2V * - 2A * - 4p s} - 

^T 1 + ir 8 = -i{12S 1 -2T 1 + 12P 1 }, (A5) 

yTi - ^T s = -i{125 8 - 2T 8 + 12P 8 } , 

\a x + \a 8 = ~{-ASi - 2Vt - 2A X + 4P} , 

y A a - ^A 8 = - i{-45 8 - 2V 8 - 2A 8 + 4P 8 } , 

\Pi + \P & = -\{Si -V 1 + \t x + A x + P} , 

yPi - \P% = -\{S 8 -V 8 + l -T 8 + A 8 + P 8 } . 

Solving these linear equations, we find that there are five independent currents. In other words, the rank of the 
10 x 10 coefficient matrix is five. Any five currents among (|A1|) are independent and can be expressed by the other 
five currents. For instance, we have the relations as 

Ss = - 7 -S 1 - 1 -V 1 - 1 -T 1+ 1 -A 1 - 1 -P 1 , 

V 8 = -2S 1 + ^V 1 +A 1 +2P 1 , 

T 8 = -6S 1 + ^Tj - 6Pi , (A6) 

A 8 = 2Si + V, + ^A x - 2P X , 

P 8 = — Si + -Vi - -T x - -A x - -Pi . 

8 2 X 2 * 4 1 2 1 6 

Finally, we establish the relations between the diquark currents and the mesonic currents. For instance, we can verify 
the relations 

^ = -1^-1^ + ^-3^-1^' 
v e = Si- ivi + \a x - Pi , 



T 3 = 


-- 35i + ir x + 3P X , 




A 3 - 


= Si + \Vi - \a x - Pi , 




P&-- 


1 1 1 1 - 

-4 5l + 4^ + 8 Tl + 4 Al - 


1 
-4 Pl 



(A7) 
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